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SSMMABX 

Is contrast to binary filters, 
continuously variable optical filters 
offer an ability to compensate for cer- 
tain inperfections in the optics of a 
hybrid correlator. Be introduce arbi- 
trary static phase errors in a nodal of 
a phase-only filtering hybrid correla- 
tor, and we siaulate a nethod of discov- 
ering a correction for then. By a re- 
cursive technique a first approximation 
to the inpulse 9 s notched filter is ad- 
justed (allowed to relax) so as to pro- 
duce successively nore localised distri- 
bution of the output in the correlation 
plane. the nethod is notivated by the 
developnent of cont inuously-variable 
phase-only spatial light nodulators, but 
it is applicable to anplitude nodulators 
and — with appropriate nodification — 
to binary nodulators as well. The tech- 
nique is robust against the fora of the 
systea* s departure froa ideal behavior. 

IBTiODUCTigil 

In a coherent optical correlator 
the iaage controlling the input spatial 
light nodulator (SLM) causes a distur- 
bance of the coherent reading wavefront. 
A lens places the optical transfora of 
the disturbance at the location of a 
second SLM, the filter plane. In clas- 
sical watched filtering, the filter has 
a nult iplicat ive effect proportional to 
the coaplex Fourier transfora of a ref- 
erence object whose correlation with the 
input pattern is desired. The coaplex 
conjugate coaprises two paraaeters, 
phase and aaplitude, but a one-paraaet er 
combination of phase and aaplitude is 
all that is usually induced with ordina- 
ry devices. The challenge is to optiaize 
a filter within the constraint of con- 
trolling only the one-paraaeter combina- 
tion of phase and aaplitude. 

Modeling of the POF (phase-only 
filter) has not yet addressed the sensi- 
tivity of the POF correlator to random 
phase deviations of the order easily 
introduced by real optics. Bsrtelt and 
Horner (1) indicate an iterative scheae 
which optimises a POF derived froa a 
classical matched filter. However, 
their technique operates entirely in 


sinulution of ideal abilities end limi- 
tations of the POF, as opposed to iter- 
ating on the basis of physical observa- 
tions. In contrast, the technique re- 
ported here is intended precisely to 
compensate for unknown, indirectly ob- 
servable departures from physical ideal- 
ity. The one-parameter form of the 
nethod is not dependent on the type of 
SLM since its relaxation is driven by 
only the input-output relationship and 
is blind to physical details of inter- 
vening interactions. 

ElACXlSdl COKSIgmXlOHS 

The pattern brought to the filter 
plane of a linear space-invariant (LSI) 
optical correlator system is not exactly 
the mathematical Fourier transfora of 
the original input image. At the filter 
plane we find the nathenatical Fourier 
transfora convolved with the diffraction 
pattern of a single pixel of the input 
spatial light nodulator, so the optical 
transfora of the LSI systea will differ 
froa 4 the mathematical Fourier transfora 
even without consideration of the non- 
ideal characteristic of the optics. In 
perforaing natched filtering or other 
correlator functions, though, we clearly 
nay need to take into account the trans- 
fora that the systea actually does on 
the input image. In addition to the 
diffraction pattern convolution, which 
does not impair the LSI nature of the 
correlator, non-ideal optical properties 
of real components and setups arise. 
Undesirable but real effects that do 
change the LSI properties of a correla- 
tor ioclude non-planarity of surfaces, 
scattering centers such as dust or in- 
clusions, and imperfect optical align- 
ment. With respect to phase, realistic- 
ally sized imperfections in a high-index 
transmitting element can cause appreci- 
able departure froa what perfect optica 
would give. A beaaspl i tting cube night 
have a tenth-wavelength non-planarity 
per surface. The number of surface 
interactions is configuration dependent, 
but three to six interactions and a 
typical refractive index give the poten- 
tial of several radians of phase distor- 
tion. For a phase-sensitive correlator 
there is clearly the possibility of 
significant unexpected distortion. 
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Filurt 1 shows ths amplitude pat- 
tsrn in ths correlation plans of an 
idealised optical correlator. A cen- 
tered inpulse is filtered with the phase 
part of its NF (notched filter). In 
Figure 2 we see the effect of adding a 
two-dineneionel seven-cycle sine-wave 
disturbance to the filtar, the distur- 
bance of the phase being only l.S ra- 
dians at its maximum. The departure 
fron the ideal notched filter is clearly 
evident. The correlation peak becones 
hardly distinguishable fron the side- 
lobes induced by the sinusoidal distur- 
bance. 



1. Correlation plane, undisturbed 
filter. 



2. Correlatioo plane, disturbed filter. 

The filtering ELM is asked to nske 
those adjustsents to a coaplex wavefront 
that will cause the wave corresponding 
to the reference isage to converge at 
the center of the correlation plane. It 
has to sccoaaodate for both the trana- 
foraing and ret raoi f oraing optica in 
order to do that. Modeling of phase- 


only notched filtering [2,3] shows that 
the POF has the possibility of producing 
very sharp correlation peaks; conconi- 
tantly one expects that phase-only cor- 
relation will be very sensitive to phase 
errors at the filter plane. 

An AOF produces broader correlations 
than the POF [3]; the hypothetical 
tenth-wave flat beansplitter affects the 
transfomed wavefront's amplitude less 
strongly than its phase; and measure- 
nents of aaplitude (intensity) are easi- 
er than neasurenents of phase. Conse- 
quently, producing a near-optinun con- 
tlnuously-vsriable filter is expected to 
be nore difficult for the POF than for 
the AOF. However, the pronise of the 
POF*s strong signal-to-noise ratio Jus- 
tifies attenpts to realise its poten- 
tial, and the easier AOF can also bene- 
fit froa the method. 

IILAXATIOH II GBA22I5X SIabch 

Suppose that the input is the de- 
sired reference iaage. Then we express- 
ly wish to have the filter that produces 
a tightly confined bright spot in the 
correlation plane. We adjust a starting 
filter so as to confine the correlation 
spot lore and lore tightly, until there 
is no adjustment of the filter that 
further confines the correlation pat- 
tern. We thus achieve the filter that is 
best matched to the reference image in 
the context of the individual correla- 
tor, including its blemishes. (If the 
correlator is not LSI, the correlation 
process is space-variant and the adap- 
tively determined filter is optimized 
only for the given position of the input 
iaage during the relaxation.) 

Central to the relaxation technique 
is a search among all possible varia- 
tions froa the current estiaate of the 
optimal filter. We need two tools. The 
first is a scalar whose value is to be 
maximized, corresponding to moving cor- 
relation plane's light into a small 
correlation spot. In practice with real 
correlators we will use the aaount of 
light falling within a snail central 
area in the correlation plane, though in 
digital aiaulation we were able to use 
sinply the aaplitude of the light at the 
central correlation spot. The second 
tool is an appropriate representation of 
the filter, s space in which to take the 
gradient of the scalar. For the space 
in which to represent the filter, we can 
use any of several candidates. In lab- 
oratory work we will use the Hadaaard 
basis, though in the simulation shown 
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here the unitary basis was used. That 
ie, we directly adjusted each pixel of 
the filter (subject to retaining POF 
characteristics), which was possible 
because of the relatively snail sise of 
the array and the 32-bit precision of 
the simulation. In the uae of another 
basis such as the Badaaard or Fourier, 
the space in which the gradient is taken 
is the one in which the coefficients of 
the filter ere the conponents. 

In general the space in which the 
gradient of the scalar is taken will 
have very high dimensionality. For an 
exhaustive search in a space of high 
dimension, the scalar must he rapidly 
calculable, and we must also consider 
the possibility of capture by false 
(local) maxima. In the relatively sim- 
ple case shown here we maximized the 
amplitude of the central correlation 
spot and, since we could evaluate the 
results directly, we did not algorith- 
mically avoid capture by false maxima. 
In a later report we will show results 
of a simulation which includes the pos- 
sibility of capture by false local maxi- 
ma. 

RESULTS 

In rigure 3 we see the results of 
spending 10 cpu hours of a TAX 11/765 
adjusting the disturbed POF that gave 
the correlation pattern of Figure 2. No 



3. Relaxed filter’s results in correla- 
tion plane. 

■ EEieri inforsatioo was given the re- 
laxation prograa regarding the fora of 
the disturbance. The prograa repeatedly 
marched in order through all the coef- 
ficients; for each coefficient the value 
of the central correlation plane ampli- 


tude was maximized. The effect of the 
adjustment order is apparent in the 
asymmetry of the correlation plane 
shape. Iteration will generally be 
required in this sort of optimisatioa 
program, but • practical implementation 
might aaek a batter order in which to 
adjust the coefficients, baaed on the 
differences between the ideal pattern 
and the current pattern. Optical corre- 
lators typically transform to amplitude 
but sense in intensity. nonetheless, 
the gradient of the intensity of the 
correlation pattern, as measured In the 
representation space of the filter, can 
be expected to have a large component in 
the direction of the inverse transform 
of the difference between the observed 
correlation intensity pattern and the 
ideal pattern. This expectation should 
lead to an optimized gradient search 
pattern. 

The non-optimal adjustment order 
wouldn’t have been so nicely visible had 
we used the Badamard or Fourier basis to 
represent the filter. It would also 
have been specious to uae the Fourier 
basis, since the disturbance has s par- 
ticularly simple representation in that 
basis and s single coefficient adjust- 
ment would have compensated for it. 

There is a significant obstacle in 
laying out an equivalent of the gradient 
search for the discrete filter. The 
obvious approach is to use the contin- 
uous representation of the filter and 
threshold it, but thresholding will 
cause an insensitivity or lack of uni- 
formity to the change in the representa- 
tion that is at the heart of the gra- 
dient search. Convergence to the opti- 
mum discrete filter is not so straight- 
forward as it is for the continuous 
filter as presented here. 

future WORK 

The work reported here was done 
with a simple low frequency form of 
phase disturbance; in the future the 
disturbance will be replaced by one with 
higher frequency (narrower autocorrela- 
tion). We will evaluate the approach to 
robustness to capture by local maxima in 
the more complicated environment, and 
the algorithm for apeading the gradient 
search will be implemented. 
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